Boundary feedback stabilization by piecewise 
constant time delay for the wave equation 

Martin Gugat 
Lehrstuhl 2 fur Angewandte Mathematik, 
Martensstr. 3, 91058 Erlangen, Germany * 

Marius Tucsnak 
Inst it ut Elie Cartan Nancy (Mathematiques) 
Universite Henri Poincar Nancy 1 
B.P. 70239, F-54506 Vandoeuvre-les-Nancy Cedex 

France ' 



Abstract 

For vibrating systems, a delay in the application of a feedback control 
can destroy the stabilizing effect of the control. In this paper we consider 
a vibrating string that is fixed at one end and stabilized with a boundary 
feedback with delay at the other end. 

We show that for certain feedback parameters the system is expo- 
nentially stable with constant delays of the form 4L/c, 8L/c, 12L/c ... 
Moreover, we show that the system is exponentially stable with piecewise 
constant delays that attain the values 4L/c and 8L/c. 
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1 Introduction 

Datko ct al. have described the following problem in the application of feedback 
laws: Some second-order vibrating systems cannot tolerate small time delays in 
their damping (see [7]). In other words: Delays can destabilize a system that is 
asymptotically stable in the absence of delays (see [5]). The problem of insta- 
bility caused by small constant delays has also been considered in [T8] , where a 
systematic frequency domain treatment of this phenomenon has been given and 
examples for the instability created by small delays have been presented. 
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In P3] a constant delay with the value 2L/c has been considered and it has 
been shown that with this delay, exponential damping is possible for feedback 
parameters with sufficiently small absolute value that have an opposite sign as 
the parameters that generate exponential damping in the case without delay. 

In this paper we show that for a constant delay that is an integer multiple 
of 4L/c, exponential damping with feedback parameters of sufficiently small 
absolute value is possible if the feedback parameters have the same sign as the 
parameters that work in the case without delay. 

Moreover, we consider piecewise constant delays with values 4L/c or 8L/c 
and show that also for delays that switch between those two values in an arbi- 
trary way, the energy decays exponentially for certain feedback parameters. 

For the problem considered in this paper some progress has been made in 
[3] for the wave equation. In [4] the related problem for the Eulcr-Bcrnoulli 
beam has been considered. In most studies of feedback stabilization of second- 
order vibrating systems, no delays are considered: In [5], a vibrating string is 
considered and a feedback law is presented for which the energy vanishes in 
finite time. In [TT] it is shown that the result from [5] is stable in the sense that 
also with moving boundaries, the energy is driven to zero in finite time. The 
problem of boundary control of the wave equation has also been studied in [20] , 
PH, [33], PH], [2], [21] and the references therein. 

This paper has the following structure: In Section[2]we define the considered 
system and in Section [3] we show that it is well-posed. 

In Scction|4]we show that the system is stable with piecewise constant delays 
that attain the values 4L/c and 8L/c. To our knowledge, this is the first example 
of a system that is stabilized with a switching delay, where the switching occurs 
between the two delay values. 

In the last section we show that our feedback law is stabilizing without delay 
and for a certain sequence of constant delays with appropriately chosen feedback 
parameters of the same sign. We show the exponential decay of the energy in 
the system. 



Let a string of length L > and the corresponding wave speed c > be given. 
Define the set O = (0, oo) x (0, L). Define the set of initial states 



Let a number t € {0, 1, 2, ...} be given. Assume that S is a piecewise constant 
function with 5{t) € [2^,4t£] for all t > 0. 
For (3/0,2/1) G B we consider the system Si : 



2 The System 



B = {(l/o,2/i) € H\0,L) x L 2 (0,L) : y (Q) = 0}. 



v(0,x) = y (x), 



(2.1) 



v t {0,x) = yi(x), x e (0, L) 
v tt (t,x) = c 2 v xx (t,x), (t,x) € O 



(2.2) 
(2.3) 
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v(t,0) = 0, t > (2.4) 

v x (t,L) = Q, *€(0,4t-) (2.5) 
c 

and 

« a (t, L) = — v t (t — 5(t), L),t> 4i-. (2.6) 
c c 

Here / is a real number. The quotient - appears in the feedback law (|2.6p in 
order to make the size of the stabilizing feedback parameters / independent of 
c. We assume that the compatibility condition yo(0) = is satisfied since it 
assures that the system has continuous states, as we show in Theorem [TJ 

3 Well-posedness of the system Si 

In this section, we study the well-posedness of system Si that is (|2.1[) - (|2.6|) . 



Theorem 1 Assume that 8 is a piecewise constant function with 6 (i) G [2 — ,4t— ] 
for all t > 0. 

Let (yo,yi) € B be given. Define the function a recursively by 

a (x) - J -|yo(-z) + yi(s)ds, x€[-L,0), ,„ ^ 

[) " I $Vo(x) + i-J X yi(s)ds, x€[0,L), (6A) 

and for k e {1, 2, 2(t - 1)} and x e [L + 2fcL, 3L + 2kL) by 

a'(x) = -a'[x - 2L) (3.2) 

and for k G {0, 1, 2, ...} ant! x € [L + AlL + 2fcL, 3L + 4ti + 2kL) by 

a'(x) = -a'(x - 2L) + fa'(x - cS(-)) - fa'ix -2L- cS(-)) (3.3) 

c c 

and the condition that a is continuous on the interval [— L, oo). Lei 

u(i,x) = a(ct + x) - a(ct- x), (t,x)ett. (3.4) 

For every finite interval I c [— Z, oo) we have a' € L 2 (I). The function v is 
continuous on fl and Vt, v x € L i 1 oc (r2). Define the family of test functions T as 

T = C 2 (S1) : There exists a set Q = [ti, t 2 ] x [a;i,X2] C SI 

such that the support of tp is contained in the interior of Q}. 

The function v satisfies the wave equation \2. S\) in the following weak sense: 

vt(t,x)(pt(t,x)d(t,x) — c 2 / v x (t,x)ip x (t,x) d(t,x) for all ip£T. (3.5) 
o Jn 

The function v satisfies \2. 1\) and 112.2]) and ^2.4^ - (27B\) . In this sense, v is the 
solution of the system Si that is 12. l\) - ![2~6\) . 
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Proof. Since y' £ L 2 (0, L), the Sobolev imbedding Theorem implies that yo 
is continuous. Moreover, y\ is in L 2 (0, L), thus a is well defined. Now we discuss 
the regularity of a. On the intervals [— L,0), [0,L) and [L + 2kL,3L + 2kL) 
(k £ {0, 1, 2, 3, ...} the function a is continuous. Due to the definition of the set 
B we have 



lim a(x) = -(l/2)j/ (0) = = (l/2)jto(0) = Km a(x), 

x— >0— x->0+ 



1 1 f L 



o 



i 



L 



yi(s) ds — a(— L) = lim a(x), 

C Jq X^L + 



a(3L + 2kL) = (f-l)a(L + 2kL)-fa(2kL-L) + C k 
= lim a(x) 

x^3L+2kL- 



hence a is continuous on the interval [-L, oo). The derivative a' in the sense of 
distributions exists on the intervals (-L, 0), (0, L), (L, 3L) and (3L + 2kL 7 5L + 
2kL) as L 2 -function. Since a is continuous, this implies that a is absolutely 
continuous on (— L, oo). Hence a' £ Lf oc (—L 7 oo). The continuity of v follows 
from the continuity of a. For t — and x £ (0, L) we have 



v(0,x) = a{x) — a(—x) — yo(x). 

For (t, x) £ almost everywhere, we have 

v t (t, x) = c[a'(x + ct) - a'{-x + ct)). (3.6) 

Thus the definition of a implies the equation Vt(0, x) = yi(x). Hence the initial 
conditions (I2.ip and (|2.2[) are valid. 

For (t, x) £ Q, almost everywhere, we have 

v x (t, x) =a'(x+ ct) + a'(-x + ct). (3.7) 



By Tonelli's Theorem (see e.g. [19]), (|3 . T[) implies v x £ L\ oc {VL) and (|3.6p implies 
vt&Ll c (Q). 
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For all <p <E T, integration by parts, p.7[) and (|3.6p yield 



v x (t,x)(p x (t,x) d{t,x) 

n 

I (p x (t, x) [a' (x + ct) + a' (— x + ct)] dt dx 

X\ J t\ 

X2 rt2 

/ ip x t{t,x)[a{x + ct) + a(—x + ct)]/cdtdx 

X\ Jti 

*2 rx 2 

/ (fitx(t, x)[a{x + ct) + a(— x + ct)]/cdx dt 

ti J xi 

t2 fX2 

\ (fit(t, x)[a'(x + ct) — a'(—x + ct)\j cdx dt 

t\ J X\ 

ipt(t,x)vt(t,x)/ c 2 d(x, t) 

u 



hence (13. 5|) holds. 

For a; = we have 0) = a(ct) — a(ct) = 0, hence at x = the boundary 
condition v(t, 0) = holds for all t > 0. 

For x = L, (|3.7[) implies for £ 6 (0, 4t-^-) the equation 

v x {t, L) = a'(L + ct) + a'{ct - L) = -a' {ct - L) + a (ct - L) = 0. 

Therefore, the boundary condition (|2.5[) holds for all t e (0,4i-^). 
For t > 4tX/c, we have 

v x (t,L) = a' (ct + L) + a' (ct - L) 

= fa'(ct- cS(t)) - fa'(ct-2L- cS(t)) 

= f[a'(L + ct- c5(t)) - a'(-L + ct - cS(t))} 

= (f/c) v t (t - cd(t), L). 



Therefore, the boundary condition (|2.6p holds for all t > 4i— . 



Remark 1 Note that our system has a continuous state. Optimal boundary 
control problems for the wave equation with countinuous states have been con- 
sidered in [§$. The proof of Theorem]^ is similar to the proof of Theorem 4-1 
in \10j - Theorem [7] is a generalization of Theorem 1 in where the case 

S(x) = 2L/c has been considered. 



6 



M. Gugat, M. Tucsnak 



3.1 Transformation of the recursion to a vector recursion 



Instead of the recursion (|3.3|) we can also use the following linear system to 
characterize the solution of Si: 



/ 


a'(x) \ 


/ 


-1 








/ 


-f \ 


/ 


a'(x-2L) \ 




a'(x-2L) 




1 
















a'\x-AL) 




a'(x-4L) 







1 













c/(x-6L) 




a'(x-6L) 










1 










c/(x-cS^)) 


V 


a'(x-c5(^) J 


V 











1 


o / 


^ a'(x 


-2L-c8(£)) J 



(3.8) 

Let B 2 be the matrix in system (|3 . 8[) . Let det(AI — B 2 ) = P/(A) denote the 



characteristic polynomial of B 2 . Then we have the equation 

p f (X) = X 5 + X 4 - Xf + f. 
If 5(x) = 4, we can write Q3.3P in the form of the linear system 



/ 




a'(x) \ 


( 


-1 


/ 


-/ 





\ 


/ 


a'(x 


-2L) \ 




a'(x 


- 2L) 




1 
















a'(x 


— AL) 




a'(x 


-4L) 







1 













a'(x 


- 6L) 




a'(x 


- 6L) 










1 










a'(x 


- 8L) 


V 


a'(x 


-8L) J 


V 











1 


o / 


V 


a'(x - 


- WL) J 



Let B\ be the matrix in system ()3.9|) . Let det(AI — B\) — df(X) denote the 
characteristic polynomial of B\. Then we have the equation 

d / (A)-A 2 (A 3 + A 2 -A/ + /). 
3.2 The characteristic Polynomial 

Let j be a natural number. For a real number / we define the polynomial 

p f (t)=t^ +1 +t^ -ft + f. 

So for / = we havepo(^) = t 2j (1 + 1) with the roots (—1) (with multiplicity 
1) and zero as the second root (with multiplicity 2j). 

Lemma 1 There exists a number 8j > 0, such that for all f G (—5j,0), all 
roots of pf have a modulus that is strictly less than one. 

For the proof of Lemma [TJ we use an intermediate result. The following 
Lemma [2] states that outside a neigbourhood of (—1), all roots of pf have a 
modulus that is strictly less than one. 

Lemma 2 Let z be a root of pf with \z — (— 1)| > i_f[i ■ Then 

\z\<l, 

that is the modulus of z is strictly less than one. 
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Proof. The equation Pf(z) = implies z 2j (z + 1) — f(z + 1) + 2/ = 0. Hence 
(z 2j —/)(« + !) = —2/ which implies the inequality 



|z 2j '-/i 



2|/l 



We have 



\z\ 2j < 



< 

Hence we have \z\ < 1 and the assertion follows. 

Note that if / < for al\t > 1 we h&vep' f (t) = t 2j - 1 ((2j+l)t+2j)-f > 

0, and since Pf(l) = 2 this implies that p/ does not have a root in [1, oo). 

Now we come to the proof of Lemma [TJ 

For t < -1, we have p'Jt) = i 2j ' _1 ((2j + l)t + 2j) - f > 0, so p f is strictly 
increasing on (— oo,— 1) and since Pf(— 1) = 2/ < this implies that p/ does 
not have a root in (— oo, —1]. Since Pf(l) = 2 > 0, there exists a real root of p/ 
in (—1,1) and all real roots of pf are contained in this interval. 

For / = 0, zo = — 1 is a single root of pj. Hence there exist numbers e > 0, 
8 > and a neighbourhood U(—l) = {z E C : \z + 1| < e} such that for all 
/ G (—5,0) there exists exactly one root of pf in U(— 1). Since the complex 
roots appear in conjugate pairs, this must be a real root, hence it is in (—1, 1). 

In particular, its absolute value is strictly less than one. 

The other 2j roots of pf are all outside of U(— 1). If |/| is sufficiently small, 
for all z ^ U(— 1) we have |z — (— 1)| > e > and thus Lemma [5] implies 

| z | < 1 which finishes the proof of Lemma [TJ 

The following Lemma implies a necessary condition that must hold if all 
roots of pf have a modulus that is strictly less than one. 



\^-f\ + \f\ 

|z + l| 1/1 
(1- I/I) + 1/1 = 1. 



Lemma 3 If f < andpf(-f) > 0, </iere exiis a rea/ rooi o/p/ wii/i absolute 
value greater than or equal to one. 

Proof. We have P/(0) = / < 0. Since p(—f) > 0, there exists a real root 
G (0,-/] with = 0. Hence ^ > 1. Let Zi,...,Z2j denote the other 

roots. Then we have \z\Zi-"Z2j\ — > L Thus the assertion follows. 



Lemma 4 If f < and A G (3 — 2^2, 3 + V2), we have 

P/((A|/|)^)>0. 
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Proof. We have 
P/((A|/|)^r) 



where 
Since 



A|/| + A^|/|w-Aw/|/|w + / 
(A-l)|/| + Aw|/|2fk+Aw|/|i^ 

l/I^TT (A- 1)|/|53TT + \27TT + \WtT\f\23TT 
\f\wk h(\f\*&), 



h(t) = A^J+rt 2 + (A - l)t + A =3+1. 



(A-1) 2 -4A = A 2 -6A + 1<0 



due to our choice of A, we have h{t) > for all t e (— oo, oo) and the assertion 
follows. 



Lemma 5 If f < 



/>-; 



1 /2j — 1 



we /iai)e 



2.7 + 1 7 2j + l 



2j V2J + 1 



4.7/ + 2 



2j 



(3.10) 



2j- 1 
2. ? + 1 



2j' 



>0, (3.11) 

(3.12) 
(3.13) 
(3.14) 

and there exist three real roots of pf, one of them in the interval (—1, — f^ipj); 
another in the interval ((— (2j |/|) =7 , — |/|2j ) and the third in ((57!/!)^, !)• 



P/(-(2j|/|)^)>0, 

p/(-|/|^)<o, 

P/((^l/l)*)<0 



/ ■ \ 2j " 

For / = -^(§+r) ,wehavep f (- 
Proof. Let i = — §fqrr - We have 



2j-l ■ 
2j+l - 



0. 



Pf (t) = t 2 \l+t) + f(l-t) 

= /( 1 + Iti) + * 2j '( 1 -^) 



2j + V 



4.7 



2.7 + 1 



2.7 + 1 



> 
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which implies p. lip . 

For A £ (l,oo), wehavep/(-(A|/|)^) = |/|[-(A + 1)A^|/|^ +A- 1] > if 

i/k^- 1 ^ 1 



A+iy A 

and p.!2p follows with the choice A = 2j. 

We have p/ (-1/1*-) = |/|[-2|/|^] < and (1XT51) follows. 

For A £ (0,1), we have P/((A|/|)^) = |/|[(A + l)AW|/|if + A - 1] < if 

i/k' 1 -^* 1 



1 + A7 A 

and (13. 14)) follows with the choice A = ^j. 

Since p/(— 1) = 2/ < and p/(0) = / < and p/(l) = 2 > the assertion 
follows. 

Lemma [3] implies that we only need to consider values of / < with 
P/(-/)<0. 

For j = 1 this yields the sharper result given in Lemma El 

Lemma 6 Let j = 1. T7ien /or a/Z / £ (1 — \/2, 0) we /iave p/(— /) < 
and Pf((~ /)/5) 1 ^ 3 ) > 0. Hence there exists a root of pf in the interval 
(— /, (— //5) 1 / 3 ). TTie otter two roote we ftat/e a modulus that is strictly less 
than one. 

Proof. We have p(-/) = -f 3 + 2f 2 + f = -/[(/-l) 2 -2]. Since / € (1-^2,0) 
we have 1 < (1 — /) 2 < 2 which implies p/(— /) < 0. Lemma S] with A = 1/5 
implies p/((— /)/5) 1 / 3 ) > 0. Hence there exists a root £* of p/ in the interval 
(-/,(-//5) 1/3 )- 

If the other roots are complex conjugate, we call them z and z and have 
|z| 2 = zz = -f/U < 1. 

Now we consider the case that the other roots are real. Note that for t > 0, 
we have p'(t) > — f > so there exists nor root that is greater than i*. On the 
other hand, for t < -1 we have p'(i) > 1 - / > 0. Since p(— 1) = 2/ < 0, this 
implies that there is no root in (—00, —1] hence also in this case the absolute 
value of all three roots is strictly less than one. Hence the assertion follows. 

For the case j ' = 2 where p/ is a polynomial of degree five we only have the 
result given in Lemma [7] 

Lemma 7 Let j = 2. Then for all f £ (-81/2500,0) we have p/(-/) < and 
the roots of pf have a modulus that is strictly less than one. 

Proof. Case 1: Suppose that pf had five real roots. Then they would all be in 
the interval (—1,1) 

Case 2: Now we consider the case that p/ has two complex conjugate roots 
z and z. Since / € (—81/(4 * 625), 0), Lemma [5] implies that we have three real 
roots h, t 2 , t 3 such that -1< h < -§ < t 2 < -|/| 1/4 < ^ |/| 1/4 < *3- 
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We have 

I i2 - / ^5 |/| 5^2^- 

Hence all roots have a modulus that is strictly less than one and the first part 
of the assertion follows with Lemma [3] 

The following Lemma gives a construction of values of / for which a pair of 
complex conjugate roots of pf for j = 2 is known. If |/| is sufficiently small and 
/ < 0, the remaining three roots are all real, so they can be easily approximated 
to arbitrary precision. 

Lemma 8 Let j — 2. Let a > be given. Define the numbers q — a+ ^ 2 ^ 2a , 
R = q + \J q 2 — 4a 3 and 

f = (8a 3 + Aa 2 )R - (1 + Aa)R 2 . 

Let b = 1 + 2a, c = lab — R and d = lac — Kb. Then we have 

p f (z) = (z 2 - 2az + R) (z 3 + bz 2 + cz + d). 

In particular, pj has the roots z\ = a + \/ R — a 2 i and Z2 = a — ^/ R — a 2 i. 

Proof. We have 

(z — z\)(z — z-i) = z 2 — 2az + R. 

Hence 

(z - z x )(z - z 2 ) (z 3 + bz 2 + cz + d) = (z 2 - 2az + R) (z 3 + bz 2 + cz + d) 
= z 5 + (b- 2a)z i + (c - 2ab + R)z 3 + (d- 2ac + Rb)z 2 + (Rc - 2ad)z + Rd 

= z 5 + z 4 + (Rc - 2ad)z + Rd. 
Using the definition of / we obtain the equation 

Rd = R[Aa 2 (l + 2a) - R(l + 4a)] = /. 
From the definition of R we have R 2 — 2qR + 4a 3 = 0. Hence 
= (1 + 2a)R 2 - 2(2a 3 + 4a 2 + a)R + 4a 3 (l + 2a). 
This is equivalent to the equation 

0= (R-2a)d + Rc. 
Hence we have — 2ad + Rc ~ —Rd = —f, thus 

(z 2 - 2az + R) (z 3 + 6z 2 + cz + d) = z 5 + z 4 - fz + f = p f (z) 
and the assertion follows. 
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4 Exponential stability of system Si with piece- 
wise constant delay 

We define the energy 

E(t) = \t (d x v(t, x)) 2 + ^ (d t v(t, x)) 2 dx (4.1) 
and the energy E\ by the equation 

E 1 {t) = Y J E(t + 2j-\. (4.2) 
j'=o ^ c ' 

Note that E(t) < Ei(t). 

To show the exponential stability of Si, we use the following result: 

Lemma 9 Let A > and the function E : [0, oo) — > [0, oo) be given. Then the 
following two statements are equivalent: 

1. E decays exponentially in the sense that there exist real numbers C\, (j, g 
(0, oo) such that 

E(t) < d E(0) exp(-/it) 

for all t G [0, oo). 

2. There exist real numbers Ci > and f G (0,1) such that the inequality 

E(t + jX) < P C 2 E(0) 

holds for all t £ [0, A) and for all j € {0, 1, 2, ...}. 

Proof. First we show that[TJ implies [5] Assume that [I] holds. Then for all 
t G [0, A) and all j £ {0, 1, 2, ...} we have the inequality 

E(t + j\) < dE(0) exp(-/i(t+jA)) 

= C\ E(0) exp(— fit) exp(—[i\j) 

< CxE(0) exp(- M Aj) 

= Ci E(0) cxp(-Xfi) j 

= PC 2 E(0) 

with C*2 = C\ and / = exp(— A/i). 

Now we show that [2] implies [1] Assume that[2j holds. For j G {0,1,2,...} 
define tj = jX. For all t G [A, oo) there exists j G {0,1,2,...} such that t G 
[tj,tj+i). Hence we can write t = t j + s, with s G [0, A). 

Define 
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Then ln(/) = — A/i. Let C\ — C2 exp(A/i). Then [21 implies the inequality 

E(t)=E(s + tj) < fC 2 E(0) 

= exp(jln(/))C 2J E(0) 

= exp(— jXfj,) C\ exp(— A/i) E(0) 

= C\ cxp(— fitj) exp(— /LtA) S(0) 

< Ci E(0) exp(-(i(tj + s)) 

= Ci £(0) exp(-/ii) 

and the assertion follows. 
Theorem 2 Lei 



_ -519801 - 761^467857 
/o 303170688 " °-° 0343 - 

Assume that the delay S is piecewise constant and that for all t > we have 
5{t) £ {AL/c, 8L/c}. 

Then there exists a neighbourhood U of fo such that for all f G U System 
Si with l — 2 is exponentially stable in the sense that the energy decays ex- 
ponentially. In fact there exists a constant Cq > that is independend of the 
initial state (yo, yi) and a constant L < 1 such that for all j G {0, 1, 2, ...} and 
for all t 6 [0, 2L/c) we have the inequality 

E 1 (t + 2j-)<L>C E l (Q). 
c 

Proof. Theorem [T] states that system Si has a solution for which we can 
compute the corresponding energy defined in (|4. 1 [) as 

E(t) = [ a'(x + ct) 2 +a'(-x + ct) 2 dx 
Jo 

L 

a'(x + ct) 2 dx. 

-L 

Let h = 2L. Let Aj, i £ {1,2,3,4,5} denote the eigenvalues of the matrix B2 
from system (|3.8I) . Assume that we have |Ab| < | A4 1 < | A3 1 < |A 2 | < |Ai|. 

Note that for a = 1/36, q = a+4 L +£ 2a3 and R = q + ^ q 2 - 4a 3 we have 
fo = 4(2a 3 + a 2 )R — (1 + 4a)R 2 . Hence for / = fo, due to Lemma [8] we have 
the eigenvalues a ± \jR—a?i . Due to Lemma [5] the other three eigenvalues are 
real and can be approximated as the roots of the polynomial of degree that is 
given in Lemma [H namely 



, 19 2 723 - V467857 3244 + 5^467857 

z A H z l H z . 

18 24624 110808 
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Define the corresponding eigenvectors 



and the matrix 



v^+Af+Af+Af+Af 



V 2 = ( Si \ S 2 S3 S4 



A? 
A? 
A, 

V 1 / 

) 



Choose the functions ci(s), C2(s), C3(s), C4(s), cs(s) such that for s G (—L,9L) 
almost everywhere we have 



/ 


a'{s- 


f 4h) \ 




/ 


ci(s) \ 




c/(s- 


f- 3/i) 






c 2 (s) 




o/(s- 








c 3 (s) 




a'(s 








C4(s) 


V 








\ 


<%(*) / 



Since the matrix is invertible and a' <E Lf oc (— L,oo) this implies ci, C2, C3, C4, 
C5 in Lf oc (—L, 9L). The functions Cj are the coefficients of the representation as 
a linear combination of the eigenvectors of the matrix B 2 ■ Then for all natural 
numbers j £ {0, 1, 2, ...} due to (|3.8|) and (|3.9|) we have the representation 



/ a'(s + 4h + jh) \ 
a'(s + ih + jh) 
a'\s + 2h + jh) 
a'(s + h + jh) 
V o/(s + jh) J 



i=l 



where 



with the matrix 



/ 7i,i+i(«) \ 
72j+i(s) 
73,i+i(s) 
74,i+i(s) 
V 75,i+i(s) I 



M{s) 



72 ,j (s) 
73, j (s) 
74 ,j (s) 
V 75, j (s) / 



M(s) 



if 
if 



6(s/c) 
6(s/c) 



8L/c, 
AL/c. 



By our construction, the matrix V 2 B2V2 — D 2 is a diagonal matrix that 
contains the numbers Ai as diagonal elements. Due to Lemma [7] this implies 
that HD2II1 < 1. (In fact, we have ||D 2 ||i < 0.994.) 

Let Hi = V 2 1 BiV 2 . Then numerical computations show that ||-Hi||i < 1. 
(In fact, we have ||iJi||i < 0.997.) 
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Define L = max{||Z? 2 ||i, ||i?i||i} < 1. Then we have the inequality 



72, j(s) 

73, j(s) 

74, j(s) 
V 75,j(s) / 



<V 



( 7i,o («) \ 
72,0 (s) 
73,o (s) 
74,0 (s) 
V 75,0 (a) 7 



= 5i J - 



(4.3) 



This implies the inequality 



/ a'(s + 4h + jh) \ 
a'(s + 3h + jh) 
a'(s + 2h + jh) 
a'(s + h + jh) 
\ a'(s+jh) J 



< II ^7ij( s ) c *( s )^ 



< 



< 



< 



i=l 



Wo 



2 2 



< ll^|| 



( 7i,j(s)ci(s) \ 

72, j(s)c 2 (s) 

73, j(s)c 3 (s) 
74 j (s)c 4 (s) 

V 75, j (a) <%(*) / 



/ 7i,j(s)ci(s) \ 
72 j (s)c 2 (s) 

73, j(s)c 3 (s) 

74, j (s)c 4 (s) 
V 75,j(s)c 5 (s) j 

1/2 



2 

1/2 



>/5||V 2 || 2 L; 



( C!(S) \ 

c 2 (s) 
c a(s) 
c 4 (s) 

V C8(«) / 

/ ci(s) \ 
c 2 (s) 

c a(s) 
c 4 (s) 

V c 5(s) y 

/ ci(a) > 
c 2 (s) 



i/2 



/ 7i,, (s) \ 

72 j (s) 

73 j (s) 

74 j (s) 

V 75 j (s) y 

/ 7i,, (*) \ 

72 ,3 (s) 

73 j (s) 

74 j (s) 

V 75 ,3 (a) y 

1/2 



c 3 (s) 
c 4 (s) 

v cb(«) y 



1/2 



2 

1/2 
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Let t £ [0, lOL/c). For the energy E\ we have the equation 

r 9L 



E^t + 2] 



L 



< 



l-L 

L 
L 



< Hv 2 \\ 2 2 Li 



y!{x + ct + jh) 2 dx 

( a'(s + ct + 4h + jh) \ 
a'(s + ct + 3h + jh) 
a'(s + ct + 2h + jh) 
a'{s + ct + h + jh) 
a'(s + ct + jh) J 

( ciOO \ 

c 2 (s) 

c 3 (s) ds 
c 4 (s) 



< 4(70^(0) 

which implies the exponential decay for f — fo due to Lemma [9l Due to 
continuity, we find a neighbourhood U of Jo such that for all f £ U we have 
||£>2(/)||i < 1 and ||i2i(/)||i < 1 and this yields the assertion. 



5 Exponential stability of system Si with con- 
stant delay 

Theorem 3 For all i £ {0, 1, 2, ...} there exists a number 8 L > such that for 
all f £ {—S L , 0) System Si with the constant delay S(t) = AlL/c is exponentially 
stable in the sense that the energy decays exponentially In fact there exists a 
constant Co > that only depends on the initial state (yotVi) an d f such that 
for all j £ {0, 1, 2, ...} and for all t £ [0, 2L/c) we have the inequality 

E(t + 2j~) < fC Q E(0). 

c 

Remark 2 Note that for the corresponding feedback law without delay 

cv x (t,L)=fv t (t,L),t>0 (5.1) 
with f = —1, the energy is controlled to zero in finite time. 

5.1 Proof of Theorem [3] 

Let i £ {0, 1,2, ...} be given. Define the characteristic polynomial Pf(t) as in 
Section [3.21 with j = i. Lemma [1] states that there exists a number S L > 0, such 
that for all / £ (— 6 L , 0), all roots of pf have a modulus that is strictly less than 
one. The proof uses the fact that from (I3.3|) we get an explicit representation 
of a' . Let zi,....,Z2 t +i denote the roots of pf. 
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Theorem [T] states that system Si has a solution for which we can compute 
the corresponding energy defined in (14.1[) as 



E(t) 



-L 



a'(x + ct) 2 dx. 



Let h = 2L. For x > (4a + 1)L equation (|3.3[) yields the equation 

a'(s) + a'{s -h)- fa'{s - 2ih) + fa'{s - (2a + l)h) = 0. (5.2) 

Using the usual method for linear difference equations we obtain an explicit 
representation of the solution a' G Lf oc (—L, oo) of (|5.2[) . Choose the functions 
ci(s), c 2 (s),...,c 2t +i(s) such that for s G (-L, (2a + 1)L) almost everywhere we 
have 



/ 


a'(s) \ 






a'(s + h) 




V 


a'{s + 2ih) j 





1 \ 



Z2L+1 



\ 



-2i+l 



J 



( 0,(8) \ 

c 2 (s) 



\ c 2t+ i(s) / 



Since the matrix is invertible and a' € Lf oc (—L, oo) this implies Cfe in L 2 oc (—L, L). 
Then for all natural numbers j G 1, 2, ... we have the representation 

a'(s+jh) = ci(s)z{ +c 2 (s)z 3 2 + ... + c 2( .+i(s) z ] 2b+l 
Let Mq = max{|2i|, \z2\, •-, \z2i+i\- Then we have the inequality 

\a'(s + jh)\ < M^\ Cl (s)\ + \c 2 (s)\ + ... + |c 2t+ i(.s)|). (5.3) 
Let t G [0, 2L/c). For the energy we obtain the inequality 



E(t + 2j~) 
c 



ct {x + jh + ct) dx 

-L 
L-ct 

i 



ct (x + ct + jh) dx 



+ [ a'(x + ct~2L + (j + l)h) 2 dx 

J L-ct 



L-ct 
L-ct 



L-ct 



Ci(x + ct)z\ + ... + c 2l+ i(x + ct)z: 



3 

2t+l 



dx 



ci(x + ct- 2L)z{ +1 + ... + c 2l+ i(x + ct- 2L)z^l x 



3+1 



dx 



< 2Ml / (| Cl ( S )| + |c 2 (.s)| + ... + |c 2t+1 (.s)|) 2 ds 



-L 



which implies the exponential decay due to Lemma |9] 
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6 Conclusion 

In this paper we have considered feedback laws that use observations from the 
past for the boundary control of the considered systems that are governed by 
the wave equation. In this way, there is enough time for the processing of the 
feedback law in practice. 

We have shown that if the feedback parameters are chosen appropriately the 
feedback laws with constant delay lead to exponential decay of the energy of 
the vibrating systems if the delay is an integer multiple of 4L/c. 

Moreover, we have shown that if the delay is piecewise constant with values 
in 4L/c, 8L/c, the system also decays exponentially if the feedback parameter 
is chosen appropriately. 
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